Vortex nucleation by collapsing bubbles in Bose-Einstein condensates 
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The nucleation of vortex rings accompanies the collapse of ultrasound bubbles in superfluids. 
Using the Gross-Pitaevskii equation for a uniform condensate we elucidate the various stages of the 
collapse of a stationary spherically symmetric bubble and establish conditions necessary for vortex 
nucleation. The minimum radius of the stationary bubble, whose collapse leads to vortex nucleation, 
was found to be 28 ± 1 healing lengths. The time after which the nucleation becomes possible is 
determined as a function of bubble's radius. We show that vortex nucleation takes place in moving 
bubbles of even smaller radius if the motion made them sufficiently oblate. 
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In this Letter we establish a new mechanism of vor- 
tex nucleation in a uniform condensate. Previously, the 
nucleation of vortices in a uniform condensate has been 
connected to critical velocities 0, 0, 0, instabilities of 
the initial states 0] or to a transfer of energy among the 
solitary waves 0- Moving positive and negative Q 
ions were shown to generate vortex rings on their surface 
where the speed of sound was exceeded. Experiments 
in superfluid helium have demonstrated long time ago 
the production of quantised vortices and turbulence |6j 
by the collapse of cavitated bubbles generated by ul- 
trasound in the megahertz frequency range. The aim of 
this Letter is to analyse theoretically for the first time 
the physics of this process in the context of the Gross- 
Pitaevskii (GP) equation. Vortex nucleation by collaps- 
ing bubbles could also be studied in the context of (non- 
uniform) atomic condensates (BEC), for which the GP 
equation provides a quantitative model, thus providing 
experimentalists with a new mechanism toproduce vor- 
tices in BEC systems, alongside rotation [8J, the decay 
of solitons and phase imprinting ^(J- Moreover, our 
work illustrates a new aspect of vortex-sound interaction 
in a Bose- Einstein condensate, a topic which is receiving 
increasing attention 

We write the GP equation in dimensionless form as 



-2i 



.dtp 



V 2 V+(l-H 2 -^(x, t))ip, 



(1) 



in dimensionless variables such that the unit of length 
corresponds to the healing length £, the speed of sound is 
c = 1/V2, and the density at infinity is = IV'ocI 2 = 1. 
To convert the dimensionless units into values applicable 
to superfluid helium-4, we take the number density as p = 
2.18 x 10 28 m~ 3 , the quantum of circulation as n = h/m — 
9.92xl0~ 8 m 2 s _:L , and the healing length as £ = 0.128nm. 
This gives a time unit 2irt; 2 / k ~ lps. Whereas for a 
sodium condensate with £ rs 0.14/im, the time unit is 



about 8ns. V(x, t) is the potential of interaction between 
a boson and a bubble. We will assume that the bubble 
acts as an infinite potential barrier to the condensate, so 
that no bosons can be found inside the bubble (ip = 0) 
before the collapse. This is achieved by setting V to be 
large inside the bubble and zero outside. 

First we consider the case of a stationary spherically 
symmetrical bubble. The spherical symmetry allows us 
to reduce the problem to dimension one, so that the equa- 
tion JU for %/} = ip(r,t) becomes 



-2tyt = V" + W/r + (1 - |Vf M 



(2) 



y 2 + z 2 . Equation is numerically 



where r 2 = x 
integrated using fourth order finite differences discretiza- 
tion in space and fourth order Runge-Kutta method in 
time. Before the collapse the field around the bubble 
of radius a is stationary, ip t = 0. The boundary con- 
ditions are ip(a,t) = stating that the bubble surface 
is an infinite potential barrier to the condensate and 
ip(oo, t) = 1. The stationary solutions for various a were 
found by the Newton-Raphson iterations. The solutions 
are ip(r) = (0, 0) if r < a and i/j(r) = (R a (r), 0) if r > a, 
with the graphs of R a {r + a) for a — 1,2,10,30 given 
on Figure 1. If the radius of the bubble, a, is suffi- 
ciently large, then we can set r = a + £ and to the 
leading order get + [1 - i?(£) 2 ]i?(0 = which 

has the solution, satisfying the boundary conditions, 
= tanh(£/V2~). The total energy of the system P3, 
£ = ±f\Vvb\ 2 dV + I J (I - \ip\ 2 ) 2 dV depends on the 
radius of the bubble, and therefore, on the form of R a : 

3 poo 

£ ™ + 2 W [R' a (r) 2 + \{l-R a {r) 2 ) 2 ]r 2 dr. (3) 



The insert on Figure 1 shows the loglog plot of the energy 
vs radius of the bubble together with the linear fit. For 
a > 20 the energy depends on the radius as £ ~ 1.65 a 2 9 . 
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FIG. 1: (colour online) The plot of the amplitude of the so- 
lution around the stationary bubbles of radii a = 1 (red), 2 
(green), 10 (blue) and 30 (black) (the smaller a corresponds 
to a steeper amplitude) . The loglog plot of the energies of the 
solutions with various a are shown on the inset together with 
the linear fit. 



FIG. 2: (colour online) The plot of the density of the conden- 
sate as a function of the distance from the centre of the cavity 
for t = (black), 10 (blue), 20 (green), 30 (red), 40 (magenta) 
after the collapse of the bubble of the radius a = 50. 
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From the energy conservation it is clear that after the 
bubble collapses and the condensate fills the cavity the 
necessary (but not sufficient) condition for vortex nuclc- 
ation is that the energy has to be greater than that of one 
vortex ring. The minimal energy of the vortex solution 
was found in to be about E ~ 55 ± 1 which corre- 
sponds to the minimum radius of a — 2.2 with £ = 55.7. 
As the condensate fills the cavity, most of the energy will 
be emitted via the sound waves, so the energy of the 
bubble has to be sufficiently greater than the energy of a 
single vortex ring to allow for such an emission. 

The time-dependent evolution of the condensate after 
the bubble collapses involves several stages. The overall 
picture is complicated by a complex interplay between 
dispersive and nonlinear effects. Dispersive effects be- 
come important on the wavelengths of order of the heal- 
ing length with the group velocity approximately given 
by d(y / k 2 /2 + fc 4 /4)/<9fc for the perturbation propagat- 
ing along the uniform state if) = 1 towards infinity and 
with the group velocity approximated by d(k 2 — l)/2dk 
for the perturbation moving along the uniform state 
ip = towards the centre of the cavity. The wavetrain 
generated by the nonlinearity is moving slower with the 
larger wavelengths than the dispersive wavetrain. Dur- 
ing the first stage dispersive and nonlinear wavetrains are 
generated near the areas of the largest curvature of the 
stationary profile. The Fourier components propagate 
at different velocities generating wave packets moving in 
opposite directions. This stage of the evolution is char- 
acterised by a flux of particles towards the centre of the 
cavity as the oscillations of the growing amplitude are be- 
ing formed on the real and imaginary parts of the wave 
function and the slope of the steep density front is get- 
ting smaller; see Figure 2. To follow the evolution of the 



particle flux we calculated the density per unit volume 
averaged over spheres of different radii, centred at the 
origin, as functions of time. Figure 3 shows these func- 
tions for the radius of the cavity a — 128 and the radii 
of the spheres over which the averaging of the density 
is performed being 4,8, 16,32, and 64. From Figure 3 it 
is clear that there is a definite moment of time t* « 97 
when the density per unit volume reaches its maximum 
at the same time (and taking the same value) for two 
smallest spheres. This moment indicates the start of a 
qualitatively new stage of the evolution in which there is 
an outward flux of particles as the condensate that over- 
filled the cavity began to expand. The density gradually 
approaches the uniform state p = 1. 

It is after the time t* that we expect the instability to 
set in and to give rise to vortex nucleation as the outward 
flow can support outward moving vortex rings. We cal- 
culated t* for various radii of the bubble and determined 
that for a > 10 the time of the start of the outflow from 
the centre of the cavity is approximated quite well by the 
linear function t* — 1.96 + 0.72 a. 

As it seen on the insets of Figure 2, some of the sur- 
faces of zero real and zero imaginary parts may be just 
a healing length apart (for instance, at t = 100 for the 
radius of the cavity a — 128, one of the surfaces of zero 
of the real part has radius 16.02 and the nearest surface 
of the zero of imaginary part has radius 17.11). Small 
radial perturbations on these zero surfaces can lead to 
their overlap which creates topological zero curves of the 
wavefunction which quickly adjust their form by emitting 
sound waves to become axisymmetrical vortex rings. 

The smaller a is, the larger the distances between the 
zero surfaces are during the condensate expansion, so it 
would require a much larger perturbation to bring these 
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FIG. 3: (colour online) The plots of the density per unit vol- 
ume as function of time for various radii of spheres over which 
the averaging is performed. The initial radius of the cav- 
ity in this case is a — 128. The radii of averaging spheres 
are b = 4,8, 16, 32, 64. The average density is calculated as 
p b = 3j^p(rydr/b s . 
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surfaces to intersect. In this case the instability mech- 
anism would be somewhat different taking longer time 
to develop. The radial 'dips' of the density of the ex- 
panding condensate, noticeable in Figure 2, are unsta- 
ble to non-spherically symmetrical perturbations, simi- 
lar to the instability of the Kadomtsev-Petviashvili 2D 
solitons in 3D [j. Depending on the energy carried by 
these 'dips,' they evolve into either vortex solutions or 
sound waves. From these considerations we expect three 
possible outcomes after bubble collapses: (1) if the ra- 
dius of the bubble is smaller then some critical radius 
a*, the density 'dips' generated by the expanding con- 
densate have rather small amplitude that decreases even 
further as they travel away from the centre quickly be- 
coming sound waves before the instability has time to 
develop; (2) after the collapse of a bubble of an inter- 
mediate size, say, of the radius a* < a < a, the waves 
of sufficiently large amplitudes are generated and the in- 
stability of these waves develops in time inversely pro- 
portional to the radius a; (3) if the radius is sufficiently 
large, a > a, the time of the first vortex nucleation is 
approximately given by the moment of the start of the 
outward flux of the particles t* approximated above. As 
a condensate continues to expand the instability mecha- 
nism described in (2) is further facilitated by the broken 
symmetry resulting from the previous nucleation events 
which leads to even more vortex rings being nucleated. 

To confirm the scenario outlined above we performed 
full three-dimensional calculations for cavities of vari- 
ous radii in a computational box of the side 200 healing 
lengths . We determined that there is a critical radius 
of the bubble for which vortex ring nucleate a* ~ 28 ± 1. 
The borderline radius between regime (2) and (3) was 
found as a ~ 45. Figure 4 shows the density isoplots 



FIG. 4: (colour online) Time snapshots of the density isoplots 
p — 0.1 after the collapse of the bubble of radius a — 50. The 
side of the computational box is 200 healing lengths and the 
distance between ticks on the side corresponds to 20 healing 
lengths. The vortex rings nucleate at about t ~ 40 (compared 
with t* ~ 38). 




at the various time snapshots after the bubble of radius 
a = 50 collapsed. Notice, that the non-symmetry of the 
field that was created at the time of nucleation (t w 40) 
continues to produce even more vortex rings as conden- 
sate expands. Each vortex ring at the moment of its birth 
has zero radius, as the surfaces of zero real and imagi- 
nary parts touch each other, and gradually evolve into a 
vortex ring of increasingly larger and radius, as clearly 
seen on Figure 4. The process in which solitary waves 
evolve into states of a higher energy was elucidated in 
0. A finite amplitude sound wave that moves behind a 
vortex ring transfers its energy to it, allowing the vor- 
tex ring to grow in size. The radius of the vortex ring 
stabilises only when it travelled sufficiently far from the 
center of the collapsing bubble, where the flow became 
almost uniform. The larger the radius of the bubble, the 
more finite amplitude sound waves will be generated at 
shorter distances, the larger the size of the final ring is 
going to be. 

So far we considered the collapse of the stationary bub- 
ble, where the vortex nucleation is connected to the in- 
stabilities developed in the spherically symmetric flow. 
There are situations when the nucleation is facilitated 
by an initial lack of the symmetry in the flow as in the 
case of a moving bubble or a bubble in the nonuniform 
(trapped) condensate. The surrounding helium exerts a 
net inward pressure across the surface, which is balanced 
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by the pressure inside the bubble. In was shown in |3( 
by asymptotic analysis of the GP equation coupled with 
the equation of the motion for the wavefunction of an 
electron that a moving bubble becomes oblate in the di- 
rection of its motion. This flattening is created by the 
difference in pressure between the poles and equator asso- 
ciated with the greater condensate velocity at the latter 
than at the former. How oblate the bubble becomes dur- 
ing its motion will depend on the velocity and pressure 
inside the bubble. The non-uniformity of the flow in the 
collapsing oblate bubble leads to vortex ring nucleation 
for bubble sizes much smaller than in the case of a sta- 
tionary spherically symmetric bubble. Figure 5 shows 
snapshots of the density plots of the cross-section of the 
collapsing bubble that prior to t — was moving with a 
constant velocity U = 0.2 and acquired an oblate form 
given by x 2 + \(y 2 + z 2 ) = 100. As the result of bubble's 
collapse four vortex rings of different radii were created. 
Three of them are moving in the same direction as the 
bubble before the collapse and one vortex ring is moving 
in the opposite direction. In the model used, the oblate- 
ness of a bubble and the velocity of its propagation, U, 
are independent parameters. A useful quantity that can 
be determined from further numerical simulations is the 
critical radius of the bubble as a function of the oblate- 
ness and U. This will be further analysed in our future 
research. 

In summary, we have suggested and studied a new 
mechanism of vortex nucleation as a result of the col- 
lapse of stationary and moving bubbles in the context of 
the GP equation. This is related to the experiments in 
superfluid helium in which the cavitated bubbles are gen- 
erated by ultrasound in the megahertz frequency range. 
Our findings suggest that a sufficiently large or deformed 
stationary bubble in a trapped condensate will produce 
vortex rings as a result of its collapse. This will also be 
a subject of our future studies. 
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